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ABSTRACT 

Under conditions in the primordial solar nebula and dense interstellar clouds, small 
grains have low relative velocities. This is the condition for efficient sticking and formation 
of fractal aggregates. A calculation of the ratio of cross-section, a , to number of primary 
particles, N, for fractal clusters yielded In a/N = 0.2635 + 0.5189N(~ 0,1748 ) . This ratio 
decreases slowly with N and approaches a constant for large N. Under the usual assumption 
of collisions producing spherical compact, uniform density aggregates, a/N varies as 
and decreases rapidly. Fractal grains are therefore much more closely coupled to the gas 
than are compact aggregates. This has a significant effect on the aerodynamic behavior of 
aggregates and consequently on their evolution and that of the nebula. 
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It has been known for some time (Whitlaw-Gray and Patterson, 1932; Fuchs, 1964) 
that aerosols are fluffy, low-density aggregates. This is also true of particulate clouds in the 
laboratory (Stephens and Russell, 1979; Samson et a/., 1987) and is illustrated in Figure 1. 
Micrometer-size or smaller grains in the atmosphere (Fuchs, 1964) or solar nebula (Volk 
et a/., 1980; Weidenschilling, 1984) have relative velocities about one ms" 1 , Therefore, 
similar structures should be found in each case. 

Recent work, both experimental (Forrest and Witten, 1979; Martin et a/., 1986) and 
theoretical (Sutherland, 1967; Witten and Sander, 1981; Kolb et al . , 1983; Meakin, 1983, 
1984a; Sander, 1984; Kolb; Samson et al ., 1987) has shown that these aggregates are frac- 
tals (Mandelbrot, 1982). Figure 2 compares a soot particle and a numerical simulation 
of cluster-cluster growth with linear trajectories. A principal characteristic of fractal ag- 
gregates is that the number of primary particles, N, within a radius r, measured from 
randomly selected particle, is given by 


N(r) = Ar°. (l) 

A is some constant dependent on the system and D is the fractal dimension which depends 
upon the accretion process. For aerosols the primary process would be cluster-cluster 
accretion (Meakin et al., 1985). This process leads to a value of D with 1.7 < D < 2.1 
(Meakin, 1984b; Julien et al., 1984). For a fractal, the density p(r) decreases with size 
according to the relation 

p(r)a r D " 3 . (2) 

For aerosol aggregates, therefore, p( r) varies approximately as r” 1 . 
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This result will have a significant effect on the aerodynamic behavior of fractal ag- 
gregates. In the usual treatments of the grain collisions where it has been assumed that 
collisions produce compact bodies with a constant density, the gas drag decreases rapidly 
with size. For fractals on the other hand, drag decreases very slowly, as we now show. 

We report here the first step towards a quantitative analysis of the aerodynamic 
behavior of fractal aggregates. A calculation of the projected cross-section as a function 
of size has been carried out. Numerical simulations were performed using the cluster- 
cluster aggregation model (Meakin et ah, 1985). A gas-kinetic velocity distribution in 
which velocities are proportional to was adopted. Nine simulations were carried out 

starting with 200,000 particles and continued until the maximum size first exceeded 10,000 
particles. 

After each cluster (containing N particles) had been generated, the projected area was 
obtained by projecting the cluster onto a plane and picking 5N points at random in a circle 
of radius R maz enclosing the projected cluster (here R max is the maximum radius of the 
cluster measured from its center of mass). If M of these points are in the region containing 
the projection, then an estimate of the projected cluster area, cx, is given by 

° = lR max ^ ( 3 ) 


For each cluster this procedure was repeated for two other projections in directions mutu- 
ally perpendicular to the first. 

The continuous line in Figure 3 shows the dependence of ln(a/N) on ln(N). Since the 
fractal dimensionality (D) is less than 2, we expect that in the asymptotic limit (N— ► oo) 
that of N should he independent of N (perhaps a could depend on some power of InfN), 
but not on a power of N). Consequently, the data shown in Figure 1 have been fitted to 
the form 

a — AN + BN^. (4) 

For clusters in the size range N — 5-100 a nonlinear least squares fit give A = 0.2635, B 
= 0.5189 and (3 = 0.8252. This curve is represented by the upper line in Figure 3 which 
has been displaced to distinguish it from the simulation data. It is apparent that this fits 
the simulation results very well for all values of N. Even for single particles (N = 1) the 
dashed curve gives a projected area of A+B or 0.782 compared to a value of n /4 or 0.785 
expected for a sphere of unit diameter. For clusters in the size range 1-2500 particles, the 
same procedure gave A = 0.2403, B = 0.5172 and j3 = 0.8465. Again, A+B = 0.768 is 
quite close to the expected values of tt/ 4. These results suggest that in the limit N — > oo 
the projected area approaches a value of about 0.24. 

For comparison, in Figure 3, we include the curve for a / N oc N -1 / 3 , appropriate for 
a compact body. This curve has been set to a = 7r/4 at N = 1 and thus matches the fractal 
curve at In N=0. Note the steep decrease in the In a/N compared to the fractal curve. 
As cr/N is the ratio of gas drag-to-particle mass, it is the dominant factor in establishing 
the gas-grain interaction. With the gas velocity field, it determines the motion of the 
aggregate. 

For fractal grains, the near constancy of this ratio means that aggregates will have 
similar aerodynamic behavior over a wide range of sizes. In particular, larger fragments 
will not settle much faster than smaller ones nor will they interact very differently with 
turbulent eddies leading to large relative velocities (Volk et al ., 1980;. Weidenschilling, 
1984). 

The optical properties of fractals also differs from those of compact particles (Berry 
and PerciVal, 1986). 
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Figure 1 - Fused Silica clutter*. Primary particle* are ~ 10 m- Because of high-formation 
temperature, particles are molten and fu»o on contact. (Courtesy of Cab-O-Sil Division, 

Cabot Corp.). 
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COMPUTER SIMULATION 


Figure 2 - Comparison of soot particle and numerical simulation at two angles of view 
showing fractal structure. (Soot-structure courtesy of G. W. Mulholland, National Bureau 
of Standards). 
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